Critical Role for 24-Hydroxylation in Homeostatic Regulation of Vitamin D Metabolism.

CONTEXT: The body has evolved homeostatic mechanisms to maintain free levels of Ca+2 and 1,25-dihydroxyvitamin D [1,25(OH)2D] within narrow physiological ranges. Clinical guidelines emphasize important contributions of PTH in maintaining this homeostasis. OBJECTIVE: To investigate mechanisms of homeostatic regulation of vitamin D (VitD) metabolism and to apply mechanistic insights to improve clinical assessment of VitD status. DESIGN: Crossover clinical trial studying participants before and after VitD3-supplementation. SETTING: Community. PARTICIPANTS: 11 otherwise healthy individuals with VitD-deficiency (25-hydroxyvitamin D [25(OH)D] ≤20 ng/mL). INTERVENTIONS: VitD3-supplements (50,000 IU once or twice a week depending on BMI, for 4-6 weeks) were administered to achieve 25(OH)D≥30 ng/mL. RESULTS: VitD3-supplementation significantly increased mean 25(OH)D by 2.7-fold and 24,25-dihydroxyvitamin D [24,25(OH)2D] by 4.3-fold. In contrast, mean levels of PTH, FGF23, and 1,25(OH)2D did not change. Mathematical modeling suggested that 24-hydroxylase activity was maximal for 25(OH)D≥50 ng/mL and achieved a minimum (∼90% suppression) with 25(OH)D<10-20 ng/mL. The 1,25(OH)2D/24,25(OH)2D ratio better predicted modeled 24-hydroxylase activity (h) (ρ=-0.85; p=0.001) compared to total plasma 25(OH)D (ρ=0.51; p=0.01) and the 24,25(OH)2D/25(OH)D ratio (ρ=0.37; p=0.3). CONCLUSIONS: Suppression of 24-hydroxylase provides a first line of defense against symptomatic VitD-deficiency by decreasing metabolic clearance of 1,25(OH)2D. The 1,25(OH)2D/24,25(OH)2D ratio provides a useful index of VitD status since it incorporates 24,25(OH)2D levels and therefore, provides insight into 24-hydroxylase activity. When VitD availability is limited, this suppresses 24-hydroxylase activity - thereby decreasing the level of 24,25(OH)2D and increasing the 1,25(OH)2D/24,25(OH)2D ratio. Thus, an increased 1,25(OH)2D/24,25(OH)2D ratio signifies triggering of homeostatic regulation, which occurs at early stages of VitD-deficiency.

Vitamin D Deficiency Increases Vulnerability to Canagliflozin-induced Adverse Effects on 1,25-Dihydroxyvitamin D and PTH.

CONTEXT: Canagliflozin has been reported to increase the risk of bone fracture-possibly mediated by decreasing 1,25-dihydroxyvitamin D (1,25(OH)2D) and increasing parathyroid hormone (PTH). OBJECTIVE: This work investigated whether baseline vitamin D (VitD) deficiency renders individuals vulnerable to this adverse effect and whether VitD3 supplementation is protective. METHODS: This community-based, outpatient study had a paired design comparing individual participants before and after VitD3 supplementation. Eleven VitD-deficient (25-hydroxyvitamin D [25(OH)D] ≤ 20 ng/mL) individuals were recruited from the Amish population in Lancaster, Pennsylvania. Participants underwent 2 canagliflozin challenge protocols (300 mg daily for 5 days): the first before and the second after VitD3 supplementation. In the VitD3 supplementation protocol, participants received VitD3 supplementation (50 000 IU once or twice a week depending on body mass index for 4-6 weeks) to achieve 25(OH)D of 30 ng/mL or greater. Two coprimary end points were identified: effects of VitD3 supplementation on canagliflozin-induced changes in 1,25(OH)2D and PTH. Secondary end points included effects of VitD3 supplementation on baseline levels of VitD metabolites and PTH. RESULTS: VitD3 supplementation increased mean 25(OH)D from 16.5 ± 1.6 to 44.3 ± 5.5 ng/mL (P = .0006) and 24,25-dihydroxyvitamin D (24,25(OH)2D) from 1.0 ± 0.1 to 4.3 ± 0.6 ng/mL (P = .0002). Mean 1,25(OH)2D and PTH were unchanged. VitD3 supplementation decreased the magnitude of canagliflozin-induced changes in 1,25(OH)2D (from -31.3%±4.7% to -9.3%±8.3%; P = .04) and PTH (from +36.2%±6.2% to +9.7%±3.7%; P = .005). CONCLUSION: VitD deficiency rendered individuals more vulnerable to adverse effects of canagliflozin on biomarkers associated with bone health. VitD3 supplementation was protective against canagliflozin's short-term adverse effects on 1,25(OH)2D and PTH.

A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies.

Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of NIPSNAP3A and an intergenic region on chromosome 1.

Type 2 Diabetes Modifies the Association of CAD Genomic Risk Variants With Subclinical Atherosclerosis.

BACKGROUND: Individuals with type 2 diabetes (T2D) have an increased risk of coronary artery disease (CAD), but questions remain about the underlying pathology. Identifying which CAD loci are modified by T2D in the development of subclinical atherosclerosis (coronary artery calcification [CAC], carotid intima-media thickness, or carotid plaque) may improve our understanding of the mechanisms leading to the increased CAD in T2D. METHODS: We compared the common and rare variant associations of known CAD loci from the literature on CAC, carotid intima-media thickness, and carotid plaque in up to 29 670 participants, including up to 24 157 normoglycemic controls and 5513 T2D cases leveraging whole-genome sequencing data from the Trans-Omics for Precision Medicine program. We included first-order T2D interaction terms in each model to determine whether CAD loci were modified by T2D. The genetic main and interaction effects were assessed using a joint test to determine whether a CAD variant, or gene-based rare variant set, was associated with the respective subclinical atherosclerosis measures and then further determined whether these loci had a significant interaction test. RESULTS: Using a Bonferroni-corrected significance threshold of P<1.6×10-4, we identified 3 genes (ATP1B1, ARVCF, and LIPG) associated with CAC and 2 genes (ABCG8 and EIF2B2) associated with carotid intima-media thickness and carotid plaque, respectively, through gene-based rare variant set analysis. Both ATP1B1 and ARVCF also had significantly different associations for CAC in T2D cases versus controls. No significant interaction tests were identified through the candidate single-variant analysis. CONCLUSIONS: These results highlight T2D as an important modifier of rare variant associations in CAD loci with CAC.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

Acute pharmacodynamic responses to sitagliptin: Drug-induced increase in early insulin secretion in oral glucose tolerance test.

Aim: DPP4 inhibitors are widely prescribed as treatments for type 2 diabetes. Because drug responses vary among individuals, we initiated investigations to identify genetic variants associated with the magnitude of drug responses. Methods: Sitagliptin (100 mg) was administered to 47 healthy volunteers. Several endpoints were measured to assess clinically relevant responses - including the effect of sitagliptin on glucose and insulin levels during an oral glucose tolerance test (OGTT). Results: This pilot study confirmed that sitagliptin (100 mg) decreased the area under the curve for glucose during an OGTT (p=0.0003). Furthermore, sitagliptin promoted insulin secretion during the early portion of the OGTT as reflected by an increase in the ratio of plasma insulin at 30 min divided by plasma insulin at 60 min (T30:T60) from 0.87+/-0.05 to 1.62+/-0.36 mU/L (p=0.04). The magnitude of sitagliptin's effect on insulin secretion (as judged by the increase in the T30:T60 ratio for insulin) was correlated with the magnitude of sitagliptin-induced increase in the area under the curve for intact plasma GLP1 levels during the first hour of the OGTT. This study confirmed previously reported sex differences in glucose and insulin levels during an OGTT. Specifically, females exhibited higher levels of glucose and insulin at the 90-180 min time points. However, we did not detect significant sex-associated differences in the magnitude of sitagliptin-induced changes in T30:T60 ratios for either glucose or insulin. Conclusions: T30:T60 ratios for insulin and glucose during an OGTT provide useful indices to assess pharmacodynamic responses to DPP4 inhibitors.

Pharmacogenetics of sodium-glucose co-transporter-2 inhibitors: Validation of a sex-agnostic pharmacodynamic biomarker.

AIM: To validate pharmacodynamic responses to sodium-glucose co-transporter-2 (SGLT2) inhibitors and test for association with genetic variants in SLC5A4, SLC5A9, and SLC2A9. METHODS: Canagliflozin (300 mg), a SGLT2 inhibitor, was administered to 30 healthy volunteers. Several endpoints were measured to assess clinically relevant responses, including drug-induced increases in urinary excretion of glucose, sodium and uric acid. RESULTS: This pilot study confirmed that canagliflozin (300 mg) triggered acute changes in mean levels of several biomarkers: fasting plasma glucose (-4.1 mg/dL; P = 6 × 10-5 ), serum creatinine (+0.05 mg/dL; P = 8 × 10-4 ) and serum uric acid (-0.90 mg/dL; P = 5 × 10-10 ). The effects of sex on glucosuria depended upon how data were normalized. Whereas males' responses were ~60% greater when data were normalized to body surface area, males and females exhibited similar responses when glucosuria was expressed as grams of urinary glucose per gram-creatinine. The magnitude of glucosuria was not significantly correlated with fasting plasma glucose, estimated glomerular filtration rate or age in those healthy individuals without diabetes with an estimated glomerular filtration rate of more than 60 mL/min/1.73m2 . CONCLUSIONS: Normalizing data relative to creatinine excretion will facilitate including data from males and females in a single analysis. Furthermore, because our ongoing pharmacogenomic study (NCT02891954) is conducted in healthy individuals, this will facilitate detection of genetic associations with limited confounding by other factors such as HbA1c and renal function.

Carriers of rare damaging CCR2 genetic variants are at lower risk of atherosclerotic disease.

Background: The CCL2/CCR2 axis governs monocyte trafficking and recruitment to atherosclerotic lesions. Human genetic analyses and population-based studies support an association between circulating CCL2 levels and atherosclerosis. Still, it remains unknown whether pharmacological targeting of CCR2, the main CCL2 receptor, would provide protection against human atherosclerotic disease. Methods: In whole-exome sequencing data from 454,775 UK Biobank participants (40-69 years), we identified predicted loss-of-function (LoF) or damaging missense (REVEL score >0.5) variants within the CCR2 gene. We prioritized variants associated with lower monocyte count (p<0.05) and tested associations with vascular risk factors and risk of atherosclerotic disease over a mean follow-up of 14 years. The results were replicated in a pooled cohort of three independent datasets (TOPMed, deCODE and Penn Medicine BioBank; total n=441,445) and the effect of the most frequent damaging variant was experimentally validated. Results: A total of 45 predicted LoF or damaging missense variants were identified in the CCR2 gene, 4 of which were also significantly associated with lower monocyte count, but not with other white blood cell counts. Heterozygous carriers of these variants were at a lower risk of a combined atherosclerosis outcome, showed a lower burden of atherosclerosis across four vascular beds, and were at a lower lifetime risk of coronary artery disease and myocardial infarction. There was no evidence of association with vascular risk factors including LDL-cholesterol, blood pressure, glycemic status, or C-reactive protein. Using a cAMP assay, we found that cells transfected with the most frequent CCR2 damaging variant (3:46358273:T:A, M249K, 547 carriers, frequency: 0.14%) show a decrease in signaling in response to CCL2. The associations of the M249K variant with myocardial infarction were consistent across cohorts (ORUKB: 0.62 95%CI: 0.39-0.96; ORexternal: 0.64 95%CI: 0.34-1.19; ORpooled: 0.64 95%CI: 0.450.90). In a phenome-wide association study, we found no evidence for higher risk of common infections or mortality among carriers of damaging CCR2 variants. Conclusions: Heterozygous carriers of damaging CCR2 variants have a lower burden of atherosclerosis and lower lifetime risk of myocardial infarction. In conjunction with previous evidence from experimental and epidemiological studies, our findings highlight the translational potential of CCR2-targeting as an atheroprotective approach.

Critical Role for 24-Hydroxylation in Homeostatic Regulation of Vitamin D Metabolism.

Context: The body has evolved homeostatic mechanisms to maintain free levels of Ca+2 and 1,25-dihydroxyvitamin D [1,25(OH)2D] within narrow physiological ranges. Clinical guidelines emphasize important contributions of PTH in maintaining this homeostasis. Objective: To investigate mechanisms of homeostatic regulation of vitamin D (VitD) metabolism and to apply mechanistic insights to improve clinical assessment of VitD status. Design: Crossover clinical trial studying participants before and after VitD3-supplementation. Setting: Community. Participants: 11 otherwise healthy individuals with VitD-deficiency (25-hydroxyvitamin D [25(OH)D] ≤20 ng/mL). Interventions: VitD3-supplements (50,000 IU once or twice a week depending on BMI, for 4-6 weeks) were administered to achieve 25(OH)D≥30 ng/mL. Results: VitD3-supplementation significantly increased mean 25(OH)D by 2.7-fold and 24,25-dihydroxyvitamin D [24,25(OH)2D] by 4.3-fold. In contrast, mean levels of PTH, FGF23, and 1,25(OH)2D did not change in response to VitD3-supplementation. Mathematical modeling suggested that 24-hydroxylase activity was maximal for 25(OH)D≥50 ng/mL and achieved a minimum (~90% suppression) with 25(OH)D<10-20 ng/mL. Suppression of 24-hydroxylase activity was triggered by mild-moderate VitD-deficiency, which sustains physiological levels of 1,25(OH)2D by suppressing metabolic clearance of 1,25(OH)2D. VitD metabolite ratios [e.g., 1,25(OH)2D/24,25(OH)2D] provide useful indices of VitD status. When VitD availability is limited, this suppresses 24-hydroxylase activity - thereby decreasing the level of 24,25(OH)2D and increasing the 1,25(OH)2D/24,25(OH)2D ratio. Thus, a high 1,25(OH)2D/24,25(OH)2D ratio signifies triggering of homeostatic regulation in response to mild-moderate VitD-deficiency. Conclusions: Suppression of 24-hydroxylase activity provides a first line of defense against symptomatic VitD-deficiency by decreasing metabolic clearance of 1,25(OH)2D. The 1,25(OH)2D/24,25(OH)2D ratio can provide a useful index of clinical VitD status.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed Whole Genome Sequencing Study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions. Large-scale whole genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess the associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with blood lipid levels (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare variant aggregate association tests using the STAAR (variant-Set Test for Association using Annotation infoRmation) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare coding variants in nearby protein coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500 kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variations and rare protein coding variations at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNA, implicating new therapeutic opportunities.

In-Depth Mass Spectrometry Analysis Reveals the Plasma Proteomic and N-Glycoproteomic Impact of an Amish-Enriched Cardioprotective Variant in B4GALT1.

B4GALT1 encodes ß-1,4-galactosyltransferase 1, an enzyme that plays a major role in glycan synthesis in the Golgi apparatus by catalyzing the addition of terminal galactose. Studies increasingly suggest that B4GALT1 may be involved in the regulation of lipid metabolism pathways. Recently, we discovered a single-site missense variant Asn352Ser (N352S) in the functional domain of B4GALT1 in an Amish population, which decreases the level of LDL-cholesterol (LDL-c) as well as the protein levels of ApoB, fibrinogen, and IgG in the blood. To systematically evaluate the effects of this missense variant on protein glycosylation, expression, and secretion, we developed a nano-LC-MS/MS-based platform combined with TMT-labeling for in-depth quantitative proteomic and glycoproteomic analyses in the plasma of individuals homozygous for the B4GALT1 missense variant N352S versus non-carriers (n = 5 per genotype). A total of 488 secreted proteins in the plasma were identified and quantified, 34 of which showed significant fold changes in protein levels between N352S homozygotes and non-carriers. We determined N-glycosylation profiles from 370 glycosylation sites in 151 glycoproteins and identified ten proteins most significantly associated with decreased galactosylation and sialyation in B4GALT1 N352S homozygotes. These results further support that B4GALT1 N352S alters the glycosylation profiles of a variety of critical target proteins, thus governing the functions of these proteins in multiple pathways, such as those involved in lipid metabolism, coagulation, and the immune response.

Acute pharmacodynamic responses to exenatide: Drug-induced increases in insulin secretion and glucose effectiveness.

AIM: Glucagon-like peptide-1 receptor agonists provide multiple benefits to patients with type 2 diabetes, including improved glycaemic control, weight loss and decreased risk of major adverse cardiovascular events. Because drug responses vary among individuals, we initiated investigations to identify genetic variants associated with the magnitude of drug responses. METHODS: Exenatide (5 µg, subcutaneously) or saline (0.2 ml, subcutaneously) was administered to 62 healthy volunteers. Frequently sampled intravenous glucose tolerance tests were conducted to assess the impact of exenatide on insulin secretion and insulin action. This pilot study was a crossover design in which participants received exenatide and saline in random order. RESULTS: Exenatide increased first phase insulin secretion 1.9-fold (p = 1.9 × 10-9 ) and accelerated the rate of glucose disappearance 2.4-fold (p = 2 × 10-10 ). Minimal model analysis showed that exenatide increased glucose effectiveness (Sg ) by 32% (p = .0008) but did not significantly affect insulin sensitivity (Si ). The exenatide-induced increase in insulin secretion made the largest contribution to interindividual variation in exenatide-induced acceleration of glucose disappearance while interindividual variation in the drug effect on Sg contributed to a lesser extent (ß = 0.58 or 0.27, respectively). CONCLUSIONS: This pilot study provides validation for the value of a frequently sampled intravenous glucose tolerance test (including minimal model analysis) to provide primary data for our ongoing pharmacogenomic study of pharmacodynamic effects of semaglutide (NCT05071898). Three endpoints provide quantitative assessments of the effects of glucagon-like peptide-1 receptor agonists on glucose metabolism: first phase insulin secretion, glucose disappearance rates and glucose effectiveness.

Vitamin D deficiency increases vulnerability to canagliflozin-induced adverse effects on 1,25-dihydroxyvitamin D and PTH.

Context. Canagliflozin has been reported to increase the risk of bone fracture - possibly mediated by decreasing 1,25-dihydroxyvitamin D [1,25(OH) 2 D] and increasing PTH. Objective. To investigate whether baseline vitamin D (VitD) deficiency renders individuals vulnerable to this adverse effect and whether VitD3 supplementation is protective. Design. This study had a paired design comparing individual participants before and after VitD3 supplementation. Setting. Community-based outpatient. Patients. 11 VitD deficient (25-hydroxyvitamin D [25(OH)D] ≤ 20 ng/mL) individuals recruited from the Amish population in Lancaster PA. Interventions. Participants underwent two canagliflozin challenge protocols (300 mg daily for five days): the first before and the second after VitD3 supplementation. In the VitD3 supplementation protocol, participants received VitD3 supplementation (50,000 IU once or twice a week depending on BMI for 4-6 weeks) to achieve 25(OH)D ≥ 30 ng/mL. Main Outcome Measures. Two co-primary endpoints were identified: effects of VitD3 supplementation on canagliflozin-induced changes in 1,25(OH) 2 D and PTH. Secondary endpoints included effects of VitD3 supplementation on baseline levels of VitD metabolites and PTH. Results. VitD3 supplementation increased mean 25(OH)D from 16.5±1.6 to 44.3±5.5 ng/mL (p=0.0006) and 24,25-dihydroxyvitamin D [24,25(OH) 2 D] from 1.0±0.1 to 4.3±0.6 ng/mL (p=0.0002). Mean 1,25(OH) 2 D and PTH were unchanged. VitD3 supplementation decreased the magnitude of canagliflozin-induced changes in 1,25(OH) 2 D (from -31.3%±4.7% to -9.3%±8.3%; p=0.04) and PTH (from +36.2%±6.2% to +9.7%±3.7%; p=0.005). Conclusions. VitD deficiency rendered individuals more vulnerable to adverse effects of canagliflozin on biomarkers associated with bone health. VitD3 supplementation was protective against canagliflozin's short-term adverse effects on 1,25(OH) 2 D and PTH.

Pharmacogenetics of SGLT2 Inhibitors: Validation of a sex-agnostic pharmacodynamic biomarker.

Aim: SGLT2 inhibitors provide multiple benefits to patients with type 2 diabetes - including improved glycemic control and decreased risks of cardiorenal disease. Because drug responses vary among individuals, we initiated investigations to identify genetic variants associated with the magnitude of drug responses. Methods: Canagliflozin (300 mg) was administered to 30 healthy volunteers. Several endpoints were measured to assess clinically relevant responses - including drug-induced increases in urinary excretion of glucose, sodium, and uric acid. Results: This pilot study confirmed that canagliflozin (300 mg) triggered acute changes in mean levels of several biomarkers: fasting plasma glucose (-4.1 mg/dL; p=6x10), serum creatinine (+0.05 mg/dL; p=8×10 -4 ), and serum uric acid (-0.90 mg/dL; p=5×10 -10 ). The effects of sex on glucosuria depended upon how data were normalized. Whereas males' responses were ∼60% greater when data were normalized to body surface area, males and females exhibited similar responses when glucosuria was expressed as grams of urinary glucose per gram-creatinine. The magnitude of glucosuria was not significantly correlated with fasting plasma glucose, estimated GFR, or age in these healthy non-diabetic individuals with estimated GFR>60 mL/min/1.73m 2 . Conclusions: Normalizing data relative to creatinine excretion will facilitate including data from males and females in a single analysis. Furthermore, because our ongoing pharmacogenomic study ( NCT02891954 ) is conducted in healthy individuals, this will facilitate detection of genetic associations with limited confounding by other factors such as age and renal function. Registration: NCT02462421 ( clinicaltrials.gov ). Funding: Research grants from the National Institute of Diabetes and Digestive and Kidney Diseases: R21DK105401, R01DK108942, T32DK098107, and P30DK072488.

Acute pharmacodynamic responses to exenatide: Drug-induced increases in insulin secretion and glucose effectiveness.

Background: GLP1R agonists provide multiple benefits to patients with type 2 diabetes - including improved glycemic control, weight loss, and decreased risk of major adverse cardiovascular events. Because drug responses vary among individuals, we initiated investigations to identify genetic variants associated with the magnitude of drug responses. Methods: Exenatide (5 µg, sc) or saline (0.2 mL, sc) was administered to 62 healthy volunteers. Frequently sampled intravenous glucose tolerance tests were conducted to assess the impact of exenatide on insulin secretion and insulin action. This pilot study was designed as a crossover study in which participants received exenatide and saline in random order. Results: Exenatide increased first phase insulin secretion 1.9-fold (p=1.9×10 -9 ) and accelerated the rate of glucose disappearance 2.4-fold (p=2×10 -10 ). Minimal model analysis demonstrated that exenatide increased glucose effectiveness (S g ) by 32% (p=0.0008) but did not significantly affect insulin sensitivity (S i ). The exenatide-induced increase in insulin secretion made the largest contribution to inter-individual variation in exenatide-induced acceleration of glucose disappearance while inter-individual variation in the drug effect on S g contributed to a lesser extent (ß=0.58 or 0.27, respectively). Conclusions: This pilot study provides validation for the value of an FSIGT (including minimal model analysis) to provide primary data for our ongoing pharmacogenomic study of pharmacodynamic effects of semaglutide ( NCT05071898 ). Three endpoints provide quantitative assessments of GLP1R agonists' effects on glucose metabolism: first phase insulin secretion, glucose disappearance rates, and glucose effectiveness. Registration: NCT02462421 (clinicaltrials.gov). Funding: American Diabetes Association (1-16-ICTS-112); National Institute of Diabetes and Digestive and Kidney Disease (R01DK130238, T32DK098107, P30DK072488).

Powerful, scalable and resource-efficient meta-analysis of rare variant associations in large whole genome sequencing studies.

Meta-analysis of whole genome sequencing/whole exome sequencing (WGS/WES) studies provides an attractive solution to the problem of collecting large sample sizes for discovering rare variants associated with complex phenotypes. Existing rare variant meta-analysis approaches are not scalable to biobank-scale WGS data. Here we present MetaSTAAR, a powerful and resource-efficient rare variant meta-analysis framework for large-scale WGS/WES studies. MetaSTAAR accounts for relatedness and population structure, can analyze both quantitative and dichotomous traits and boosts the power of rare variant tests by incorporating multiple variant functional annotations. Through meta-analysis of four lipid traits in 30,138 ancestrally diverse samples from 14 studies of the Trans Omics for Precision Medicine (TOPMed) Program, we show that MetaSTAAR performs rare variant meta-analysis at scale and produces results comparable to using pooled data. Additionally, we identified several conditionally significant rare variant associations with lipid traits. We further demonstrate that MetaSTAAR is scalable to biobank-scale cohorts through meta-analysis of TOPMed WGS data and UK Biobank WES data of ~200,000 samples.

Whole genome sequence analysis of apparent treatment resistant hypertension status in participants from the Trans-Omics for Precision Medicine program.

Introduction: Apparent treatment-resistant hypertension (aTRH) is characterized by the use of four or more antihypertensive (AHT) classes to achieve blood pressure (BP) control. In the current study, we conducted single-variant and gene-based analyses of aTRH among individuals from 12 Trans-Omics for Precision Medicine cohorts with whole-genome sequencing data. Methods: Cases were defined as individuals treated for hypertension (HTN) taking three different AHT classes, with average systolic BP ≥ 140 or diastolic BP ≥ 90 mmHg, or four or more medications regardless of BP (n = 1,705). A normotensive control group was defined as individuals with BP < 140/90 mmHg (n = 22,079), not on AHT medication. A second control group comprised individuals who were treatment responsive on one AHT medication with BP < 140/ 90 mmHg (n = 5,424). Logistic regression with kinship adjustment using the Scalable and Accurate Implementation of Generalized mixed models (SAIGE) was performed, adjusting for age, sex, and genetic ancestry. We assessed variants using SKAT-O in rare-variant analyses. Single-variant and gene-based tests were conducted in a pooled multi-ethnicity stratum, as well as self-reported ethnic/racial strata (European and African American). Results: One variant in the known HTN locus, KCNK3, was a top finding in the multi-ethnic analysis (p = 8.23E-07) for the normotensive control group [rs12476527, odds ratio (95% confidence interval) = 0.80 (0.74-0.88)]. This variant was replicated in the Vanderbilt University Medical Center's DNA repository data. Aggregate gene-based signals included the genes AGTPBP, MYL4, PDCD4, BBS9, ERG, and IER3. Discussion: Additional work validating these loci in larger, more diverse populations, is warranted to determine whether these regions influence the pathobiology of aTRH.

A framework for detecting noncoding rare-variant associations of large-scale whole-genome sequencing studies.

Large-scale whole-genome sequencing studies have enabled analysis of noncoding rare-variant (RV) associations with complex human diseases and traits. Variant-set analysis is a powerful approach to study RV association. However, existing methods have limited ability in analyzing the noncoding genome. We propose a computationally efficient and robust noncoding RV association detection framework, STAARpipeline, to automatically annotate a whole-genome sequencing study and perform flexible noncoding RV association analysis, including gene-centric analysis and fixed window-based and dynamic window-based non-gene-centric analysis by incorporating variant functional annotations. In gene-centric analysis, STAARpipeline uses STAAR to group noncoding variants based on functional categories of genes and incorporate multiple functional annotations. In non-gene-centric analysis, STAARpipeline uses SCANG-STAAR to incorporate dynamic window sizes and multiple functional annotations. We apply STAARpipeline to identify noncoding RV sets associated with four lipid traits in 21,015 discovery samples from the Trans-Omics for Precision Medicine (TOPMed) program and replicate several of them in an additional 9,123 TOPMed samples. We also analyze five non-lipid TOPMed traits.

Rare Variants in Genes Encoding Subunits of the Epithelial Na+ Channel Are Associated With Blood Pressure and Kidney Function.

BACKGROUND: The epithelial Na+ channel (ENaC) is intrinsically linked to fluid volume homeostasis and blood pressure. Specific rare mutations in SCNN1A, SCNN1B, and SCNN1G, genes encoding the α, ß, and γ subunits of ENaC, respectively, are associated with extreme blood pressure phenotypes. No associations between blood pressure and SCNN1D, which encodes the δ subunit of ENaC, have been reported. A small number of sequence variants in ENaC subunits have been reported to affect functional transport in vitro or blood pressure. The effects of the vast majority of rare and low-frequency ENaC variants on blood pressure are not known. METHODS: We explored the association of low frequency and rare variants in the genes encoding ENaC subunits, with systolic blood pressure, diastolic blood pressure, mean arterial pressure, and pulse pressure. Using whole-genome sequencing data from 14 studies participating in the Trans-Omics in Precision Medicine Whole-Genome Sequencing Program, and sequence kernel association tests. RESULTS: We found that variants in SCNN1A and SCNN1B were associated with diastolic blood pressure and mean arterial pressure (P<0.00625). Although SCNN1D is poorly expressed in human kidney tissue, SCNN1D variants were associated with systolic blood pressure, diastolic blood pressure, mean arterial pressure, and pulse pressure (P<0.00625). ENaC variants in 2 of the 4 subunits (SCNN1B and SCNN1D) were also associated with estimated glomerular filtration rate (P<0.00625), but not with stroke. CONCLUSIONS: Our results suggest that variants in extrarenal ENaCs, in addition to ENaCs expressed in kidneys, influence blood pressure and kidney function.

Whole genome sequence analysis of blood lipid levels in >66,000 individuals.

Blood lipids are heritable modifiable causal factors for coronary artery disease. Despite well-described monogenic and polygenic bases of dyslipidemia, limitations remain in discovery of lipid-associated alleles using whole genome sequencing (WGS), partly due to limited sample sizes, ancestral diversity, and interpretation of clinical significance. Among 66,329 ancestrally diverse (56% non-European) participants, we associate 428M variants from deep-coverage WGS with lipid levels; ~400M variants were not assessed in prior lipids genetic analyses. We find multiple lipid-related genes strongly associated with blood lipids through analysis of common and rare coding variants. We discover several associated rare non-coding variants, largely at Mendelian lipid genes. Notably, we observe rare LDLR intronic variants associated with markedly increased LDL-C, similar to rare LDLR exonic variants. In conclusion, we conducted a systematic whole genome scan for blood lipids expanding the alleles linked to lipids for multiple ancestries and characterize a clinically-relevant rare non-coding variant model for lipids.